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Final-state Polarization in Decays 
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Certain Eg — > V1V2 decays {Vi is a vector meson) can be related by flavor SU(3) symmetry to 
corresponding V3V4 decays. In this paper, we show that the flnal-state polarization can be 

predicted in the B^ decay, assuming polarization measurements of the 5° decay. This can be done 
within the scenario of penguin annihilation (PA), which has been suggested as an explanation of the 
unexpectedly large transverse polarization in B ^ (j)K* . PA is used to estimate the breaking of flavor 
SU(3) symmetry in pairs of decays. Two of these for which PA makes a reasonably precise prediction 
of the size of SU(3) breaking are (B^ ^ <^(^, ^ 0^"°*) and (B° ^ (l)K°*,Bl K°*K'"'). The 
polarization measurement in the -B^ decay can be used to predict the transverse polarization in the 
Bs decay, and will allow a testing of PA. 



PACS numbers: 13.25.Hw, 13.88.+0, ll.30.Er 
INTRODUCTION 

We consider B — > V1V2 decays {Vi is a vector meson). 
Since the final-state particles are vector mesons, when the 
spin of these particles is taken into account, this decay 
is in fact three separate decays, one for each polariza- 
tion (one longitudinal, two transverse). Naively, within 
the standard model (SM), the transverse amplitudes are 
suppressed by a factor of size niv/mg {V is one of the vec- 
tor mesons) with respect to the longitudinal amplitude. 
Then one expects the fraction of transverse decays, /t, 
to be much less than the fraction of longitudinal decays. 

However, it was observed that these two fractions are 
roughly equal in the decay B — > (j)K*: JtI Jl — 1 
A similar effect was later seen in i? ^ pK* decays 

If one goes beyond the naive SM, there are two expla- 
nations [5|] which account for this "polarization puzzle." 
The first is penguin annihilation (PA) Q. i3 ^ (jjK* re- 
ceives penguin contributions, bOsqOq, where q — u,d 
(O are Lorentz structures, and color indices are sup- 
pressed). With a Fierz transformation, these operators 
can be written as bO'qqO's. A gluon can now be emit- 
ted from one of the quarks in the operators, and can then 
produce a pair of s, s quarks. These then combine with 
the s, q quarks to form the final states (pK*^ {q = u) 
or (j)K*'^ {q = d). These are annihilation contributions. 
Normally such terms are expected to be small as they 
are higher order in the l/nib expansion, and thus ig- 



nored. However, within QCD factorization (QCDf) 0], 
it is plausible that the coefficients of these terms are large 
0]. (Within perturbative QCD ^, the penguin annihila- 
tion is calculable and can be large, though it is not large 
enough to explain the polarization data in _B ^ ipK* [9].) 

In QCDf, due to the appearance of endpoint diver- 
gences, PA is not calculable, but is modeled f^. These 
divergences are regulated with a cut-off, introducing sev- 
eral arbitrary parameters. There is therefore an enor- 
mous uncertainty in the size of the PA amplitude as one 
varies these unknown parameters within certain chosen 
limits 

It is also possible within QCDf that the transverse am- 
plitudes receive significant contributions from perturba- 
tive rescattering from charm intermediate states. How- 
ever, the transverse amplitudes could be purely domi- 
nated by PA. In this paper we explore the consequences 
of the scenario in which PA contributions are large and 
dominant to see what type of testable predictions result. 

The second SM explanation is rescattering 
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The idea is that nonperturbative rescattering effects in- 
volving charm intermediate states, generated by the op- 
erator bO'ccO's, can produce large transverse polariza- 
tion in _B — !■ 4>K* . A particular realization of this sce- 
nario is the following Consider the decay i?+ 
]j*+f)*o Since the final-state vector mesons are heavy, 
the transverse polarization can be large. The state 
]j*+]j*o (ng^j^ j^Q.^y rescatter to (pK*~^. If the transverse 
polarization T is not reduced in the scattering process, 
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this will lead to (j)K*~^ with large fr/ II- (A simi- 

lar rescattering effect can take place for (j)K*^.) 

It is important to test these explanations in order to 
determine whether new physics is or is not present. The 
polarization puzzle has been mainly seen in & ^ s tran- 
sitions. However, if PA or rescattering is the true expla- 
nation, one also^xpects to observe large fr/ Jl in b —^ d 
decays. In Ref. [51] such decays were discussed, and it was 
observed that the most promising transitions were those 
which are dominated by penguin amplitudes. The b ~^ s 
and b d penguin decays are 



b- 


-> sss and b - 


-> sdd 




^ 4>K°* 


















B+ - 


(t>K+*,p^K' 


b- 


-> dss and b - 


-f ddd 






















B+ - 





(1) 



(Decays which also receive tree contributions are not in- 
cluded in the current analysis.) 

Now, all of these decays are the same under flavor 
SU(3), which treats d, s and u quarks as equal. The idea 
is that, given a measurement of the polarization in one 
decay, one can predict the polarization in another decay 
using PA or rescattering. However, in relating the two 
decays, the effect of SU(3) breaking must be included. 
We can relate the transverse amplitudes of SU(3)-related 
decays in the scenario in which PA dominates these am- 
plitudes. On the other hand, this relation is unknown in 
rescattering, which involves long-distance contributions. 
For this reason, in this paper we consider only PA. 

We note that the transverse and helicity ampli- 

tudes {A±) are related by_yljj_j^ =_{A+ ± A_) /\/2. How- 
ever, A_ for B decays (^4. for B decays) has an extra 
spin-flip 0(l/mb) suppression with respect to (^-)- 
Consequently, we neglect A_ (A+) and henceforth define 
= + |A^|2 « (i.e. Ar includes both 

transverse amplitudes). 

First, consider the pair of decays i?" — 
(j>K'^* . The main PA contributions to the transverse am- 
plitudes are (the penguin-annihilation term arises only 
from penguin operators with an internal t quark) 



and B^ 



= \VtlVts\fsof, 2(6<' 



for'), 



(2) 



where 63^^ and 64^^^^-' are the QCDf terms correspond- 
ing to PA Throughout the paper, we have dropped 
the overall factor of Gf/V2. (Absolute values are taken 
for the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments because we are not interested in CP- violating ob- 
servables here, but rather the rate.) 

As noted earlier, QCDf also contains the (perturbative 
rescattering) term 04 which can contribute here. It is 



the coefficient of the {V — A) ^ {V — A) piece in the 
operator product expansion. However, the effect of 04 
could be small in the transverse amplitude because of 
spin flips. Consequently, it is reasonable to ignore 04 in 
pure penguin decays, and we do so here. 

The annihilation coefficient 63 is given by [13] 

[(C5 + N,Ce) Ai + C5AI + C3 A\] , (3) 
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N? 



where Nc is the number of colors, the Ci are Wilson co- 
efficients, and the incalculable infrared divergences are 
found in the A]!'^ (due to the endpoint singularity of the 
final-state distribution amplitudes). The superscripts 'i' 
and '/' refer to giuon emission from the initial- and final- 
state quarks, respectively. The subscript '1' refers to the 
Dirac structure {V — A) (g) {V — A), while '3' refers to the 
Dirac structure (-2)(S' - P) (g) {S + P). 
The annihilation coefhcient 64 is 



C 



4 Ar2 



l(V-iV-2) _ 



^{C,A\ + C,Al) 



(4) 



Here, A2 « A\ has a suppression factor of 0(l/m^) com- 
pared to the A3-' [T~ 



Thus, 6^^^' sind the third term 



can be neglected in Eq. 



Now, if we assume that the term containing A^ in 
^^^1^2) dominates over the others (as we will see below, 
the error on this approximation is at the level of only a 
few percent), we find that 



AM 



has the following integral form flo| 



where A^'^^^''^' 



(5) 



(6) 



dxdyi 



2mi 



V2 



T 



(t>al{x)(l)2{y) , 2m2 ^ 6i{x)(t)b2{y)\ 



mi 



x'^y 
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Here cf)^ {4'a,b) is the twist-2 (twist-3) light-cone dis- 
tribution amplitude, x (y) stands for the momentum 
fraction carried by the quark in Vi (V2), and — 
{2mv/mb)f^/fv fv is defined as [13] 



(7) 



It is useful to make a comment concerning the appro- 
priate interpretation of Eq. ©: if a K* is one of the 
final particles, the argument of the distribution ampli- 
tudes (DAs) is the momentum fraction corresponding to 
the s quark, whereas for a K* meson, it is the momentum 
fraction of the d quark. That is, the DAs are defined as 
4>K' (z) — 4>K' {z), where z =1 — z. 

In order to estimate the ratio of the transverse ampli- 
tudes in Eq. ([5]), we need to know the amount of SU(3) 
breaking in the ratio of A^***' and A^'*"^"*'. To do this, 
further assumptions are necessary. In what follows, we 
adopt the same assumptions as those used by the authors 
of Ref. fioj to carry out their study: 
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1. the asymptotic form of the hght-cone distribution 
(LCD) ampUtudes, 

2. a universal parametrization of the end-point sin- 
gularities, i.e. independent of any particular decay 
mode, 

3. a modeling of the singularities. 

It is the first point about which there has been some 
debate. Certain references have calculated higher-order 
moments in the LCDs, suggesting that such "non- 
asymptotic LCDs" are important for some light mesons 
14| . If so, then SU(3) breaking in these non-asymptotic 



pieces will also contribute to the ratio in Eq. Unfor- 
tunately, this SU(3) breaking is not calculable, in which 
case our analysis below will not hold. Equally unfortu- 
nately, it is very difficult for experiment to determine 
which type of LCD is present ^1^. Thus, the reader 
should be aware that our predictions are not only a test 
of PA dominance, but also of asymptotic LCDs. 

Now, the DAs for the final states are universal in the 
asymptotic limit, i.e. (j>'^{x) = 6xx and (t)a{x) — (f>b{x) = 
3x^. Thus, in this approximation, we find that Ag***' 
and Ag'*^ ' have exactly the same dependence on x and 
y [Eq. ([6])]. In this case, the ratio between the transverse 
amplitudes becomes 



fB° 



2rf 



where /bo//bO = 1.22 ± 0.03 [161, U = 221 ± 3 MeV 
[11, andV^o- = 218 ± 4 MeV III]. The values of 
in are theoretically estimated. We have [12] fj- = 
= 175 ± 25 MeV. Since the decay constants and the 
meson masses are known, the SU(3) breaking in Eq. ([5]) 
is well-controlled for this pair. 

What this says is that the transverse polarization am- 
plitude in (t)(j) is predicted by PA to be related to 
that in (j)K"* through Eq. Thus, once one 
makes the polarization measurement in the decay, 
one can test PA by making the equivalent measurement 
in the B^ decay. 

The key ingredient in the above analysis is to take two 
decays in which the final states have the same dependence 
on the momentum fractions x and y. However, although 
the pair of decays considered above is the most promis- 
ing for the analysis, it is not unique. In fact, all decays 
in a special class have the same dependence on x and y. 
This class contains b{B'^, B^) — > sdd (the parenthesis in- 
dicates that this transition includes only 5° or 5+, and 
not , decays) and b sss penguin decays. In addition, 
only decays to ground-state spin-1 mesons are included 
(excited mesons have different DAs in general). Thus, 
excluding those decays which also receive tree contribu- 
tions, these correspond to B+ p+K°*, B^ (f>K°*, 



_B+ 4)K** , and (fxp. The important point here 

is that the dependences on the momentum fractions in 
are the same for every decay belonging to this class. 
Therefore, in the comparison of any two of these decays, 
the integrals containing the singularities cancel in the 
ratio of the transverse amplitudes, and this even before 
using a cutoff to regulate the end-point divergences. 

Of course, the size of the SU(3) breaking will depend on 
the pairs of decays considered [see Eq. ([8])]. For example, 
we expect the pair B^ <j)(j) and B* (j)K'^* to be as 
good as B° ct)(j} and ^ (j}K°*. However, the SU(3) 
breaking in _B+ p^K^* and i?+ — > cf)K^* also turns 
out to be about the same size. In all cases, the SU(3) 
breaking can be worked out as we have done above. 

Another class in which the final states have the same 
dependence on x and y is given by the decays b{B^) 
ddd and b{Bl) —> sdd. Even if one of the most promising 
B^ decay modes to be measured in the near future, 
K°*K°*, belongs to this class, we cannot use it to make 
predictions because its decay-class partner, B° — + p°K*°, 
typically receives tree contributions. Therefore, this class 
is not very useful since it only contains one pure penguin 
decay. 

The third class is defined by the transition 6 — > dss. 
This includes the decays B'^ cj)K^* and 5° K°*K°*, 
and this pair is particularly promising. The prediction for 
the ratio of the transverse amplitudes is 



AriBO^^K^*) 



fB<>h 



(9) 



The pair 5° (j)K"* and B+ -> K+*K"* can be treated 
similarly. 

An important consequence of the above discussion is 
that it is not possible to get cancellations in the ratios 
through the comparison of decays belonging to different 
classes. For example, we will not obtain a clean result 
by comparing B^ 



K°*K°* and B° 



K°*K°\ even 



though the final states are identical in the two decays. In 
addition, we have presented the list of decays in Eq. ([T]) 
in terms of 6 — > d or 6 — > s transitions, motivated by 
the apparent predominant role played by the b s tran- 
sitions in the polarization puzzle. However, as shown 
above, this is not the most natural way to classify the 
decays in order to achieve the most reliable predictions 
for penguin annihilation within QCDf. 

We now turn to the estimation of errors in Eqs. ([8]) 
and ([9]) due to the single inclusion of A^ in the trans- 
verse amplitudes. We study the relative magnitude of 
the A\ g terms in 63 [Eq. ([3])], as well as the relevance of 
the annihilation coefficient 64 [Eq. (|4])] in the cases it is 
appropriate. To do this, we follow several applications of 
QCDf in which the incalculable infrared divergences in 
the A^f^-^ are isolated with a cutoff and the penguin an- 
nihilation amplitude is modeled by introducing unknown 
parameters. 
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In passing, we note the following. Previously, we men- 
tioned that we use the same assumptions as those in 
Ref. Although we work within the same restricted 

theoretical framework as this reference, and although it 
is true that the size of the error in our predictions could 
be affected by large uncertainties related to the choice 
of this particular scenario, we emphasize that our results 
go beyond the analysis made in Ref. [13] ■ There, due to 
the parametrization of the infinities, the uncertainties in 
the individual transverse amplitudes turn out to be at 
the level of one hundred percent for most of the decay 
channels [l^ . Instead, we show here that it is possible to 
obtain more accurate predictions with the same theoreti- 
cal inputs used in the treatment of the divergent integrals 
when specific decays are compared. 

To illustrate the procedure used in the estimation of 
errors, we focus on the first pair, (fxj) and Sj] 

(j)K^* (a similar reasoning holds for <j)K^* and 

K°*K°*, as well as several other decay pairs). 
We first write the ^{'g after the parametrization of the 
infrared divergences is applied 



oeps 

M 



mi 



^TO2 * mi ^ 



2 2 

) (X'i - 2X^ + 2) 

4 « 187ra, (— r^^ + —vY^) (2X1 - 5X^ + 3), (10) 
TO2 nil 

We see that A\ 2 and are symmetric in the interchange 
of Vi V2, while A\ is antisymmetric. X^ and X^ 
contain the same input parameters, but have different 
end-point-divergence behavior: 

AT ' 



X^ = il + p^e'^^)ln- 



X, = (l + p.e^*-)^ 
Ah 



(11) 



Here, Ah is an input parameter (A/j = 0.5 GeV [1(|), and 
(j}^ is an arbitrary phase. We have taken Ah, (/ja, and 
to be the same for every decay mode. 

To study the relative significance of 64 and the ne- 
glected terms in 63, we evaluate the following ratios 
[Eqs. © and dl]: 



„(^'l^2) 



'bi 



JV1V2) 



R 



C4 


+ Ce 


4 '(^1^2) 

^1,2 


C5- 
















j^HViV2) 




C5 




C5- 


f iVcCe 


Af(<l>S<°') 

^3 



(12) 



where V1V2 = #, (l)K^* (note that Rl^^' is zero). 

Although the values of these ratios are quite uncertain, 
in large part due to the (arbitrary) value of Ah, in vir- 



< 



l'l-"'f)3 I ^ 0(10" ). One can see this as fol- 
lows. First, the ratios of the relevant Wilson coefficients 
in Eq. (fT^ at n — mb/2 are as follows 0: 



C5 + NcCe 



C4 + Cg 

C5 + N,Ce 

« -0.11 



0.63 

C5 



C5 + NcCe 



-0.05 ,(13) 



Second, we have (Ar^^^'/A^'^^''^') 



|(a-6)/(a 



O {11111712/112^), and 
6)1 « 0.07 (where a 



denotes {rrij^o, /m^)r^ and b is given by K ^ (f>). 

We evaluate the three ratios in Eq. by considering 
many different values in the ranges < < 2 and < 
0, < 27r. We find that \r™\, |r<;-^^^'|, |i?<f"*>| « 1 
always, except for a singular behavior at = 0, 27r. The 
largest contribution to the error arises from f^^^^^' but it 
remains at the level of a few percent within the scanned 
region of the parameter space. Thus, we have covered a 
wide set of models of the infrared singularities. The point 
here is that, although the precise values of the ratios are 
very uncertain, they are always small. 

We therefore conclude that the PA dominance hypoth- 
esis leads to a clean prediction for the ratio of trans- 
verse amplitudes in the pair — > (j)(f> and B^ — > (j)K'^* 
[Eq. ([8])], though this result is a direct consequence of the 
particular modeling of the suppressed terms (e.g. asymp- 
totic LCDs). We have also analyzed the pair i?" (j)K^* 
as well as the pairs containing the 



and B^ K°*K° 
charged modes and (B"*" 



-K"*, B^ 



tually all cases it is found that 



JVlV2)\ < 



0(10"^) and 



der the same set of assumptions, and we have obtained 
similar conclusions. 

We have therefore seen that there are a number of de- 
cay pairs within a given class whose transverse polariza- 
tions are related. Some of these pairs involve a B'^ decay. 
Now, the B-factories BaBar and Belle have made many 
measurements of B^ and _B+ mesons. But it is only rela- 
tively recently, at hadron colliders, that 5° mesons have 
started to be studied. This will increase when the LHCb 
turns on. It should be possible to make measurements 
of the transverse polarization in some i?^ decays in the 
near future, and to test the PA/QCDf hypothesis. 

Above, we have presented the ratio of ^^'s for two 
pairs of decays. However, it is perhaps better to present 
a ratio of /^'s since this is what will actually be measured. 
At and frp are related by including information about the 
branching ratio (BR): = |yl^|2/(r BR/ PS), where F 
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is the total width and PS is the phase space. We find 



MB", ^ ct^K"*) 

3.22 ±0.72 

MB", ^ K"*K"*) ^ 
MB"s - ^K"*) 

0.62 ±0.12 



0*\ 



BR{B", - 

BRiB", 4>K"*) 
BR{Bl ^ K"*K"* 



.(14) 



The numbers have been obtained by taking values of 
masses and lifetimes from the Particle Data Group with- 
out errors [13], along with the theoretical estimates of 
the decay constants given above. The predictions given 
in Eq. ([14]) will yield a test of PA. If there are discrepan- 
cies in the measurements, this may indicate the presence 
of new physics, in & — > s and/or b d transitions. 

The decays B^ k"*K"* [H and ^ (j)K"* 
have both been measured, so that this information can 
be included in Eq. ifM)) : 



Bd 



0* 



B2 



K"*K' 



,PK"* 



h = o.8it°:l° ± 0.06 , 

BR = (OMtoAl) X 10"^ 
A = 0.49 ± 0.04 , 
BR= (9.5 ±0.8) X 10"^ 



(15) 



fj, is defined as = 1 — /l . 

To summarize, a large /t/ II has been observed in sev- 
eral b s B V1V2 decays. There are two expla- 
nations of this measurement within the standard model 
- penguin annihilation (PA) and rescattering. Now, one 
logically also expects to see a large Jt/Jl in certain b ^ d 
decays. The most promising decays are those dominated 
by penguin amplitudes, and there are quite a few 6 — > d 
and & — > s penguin decays. All of these are equal un- 
der flavor SU(3) symmetry. Given the measurement of 
fr/ fh in one decay, if one wishes to predict fr/ /l in an- 
other decay, it is necessary to take SU(3) breaking into 
account. However, it is only within a specific scenario of 
PA dominance for the transverse amplitudes, and for spe- 
cial classes of decay pairs, that this SU(3) breaking can 
be estimated. We therefore assume that it is PA alone 
which is the source of the large transverse polarization 
and explore its consequences. 

We find that there are several decay pairs for which 
PA makes a reasonably precise estimate of the SU(3) 
breaking (assuming asymptotic LCDs). Thus, given the 
measurement of Jt/Jl in one decay, PA makes a pre- 
diction for the transverse polarization in the second de- 
cay. In this paper we have concentrated on two decay 
pairs that involve B^ mesons: (B^ 0^,-8^ — > <pK"*) 
and (B° ^ (j)k"*,B", ^ K"*K"*). The polarization 
measurement in the B^ decay allows one to predict the 
transverse polarization in the B" decay. This will permit 



the explicit testing of PA, probably in the near future at 
the LHCb. 
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